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I. INTRODUCTION 



The critical behaviour of weakly-disordered quenched systems undergoing continuous 
phase transitions is of great interest. The study of critical properties of random spin systems 
in which the local energy density couples to quenched disorder has a long history going back 
to the classical papers by A. B. Harris and T. C. Lubensky and D. E. Khmelnitskii 
0]. They initiated a considerable progress in studying disordered systems by applying the 
conventional field-theoretical renormalization-group (RG) approach based on the standard 
scalar Lp^ theory with n-component order parameter in (4 — e) space dimensions. At present 
it is commonly believed that the RG approach provides a thorough understanding how weak 
disorder affects thermodynamic properties of random systems in a close vicinity of the Curie 
point. 

According to the Harris criterion, critical exponents of weakly disordered Ising model 
should differ from those for the pure system [Q. The first regular method for calculating 
their values, famous ^e-expansion, was invented more than 20 years ago but turned 

out to be numerically ineffective, at least in lower orders in Fair numerical estimates 
for critical exponents of the random Ising model (RIM) were obtained in the framework of 
the renormalization-group approach in three dimensions from two-loop [Q, three-loop 
and four-loop RG expansions. The RG method at fixed dimensions proved also to 
be efficient when used to calculate critical exponents of the two-dimensional RIM and the 
marginal value Uc of the order parameter dimensionality n for the cubic model P,p|, p!0| , p 



as 



is known, ric separates the region where the system becomes effectively isotropic approaching 
the critical point {n < ric) from that of the essentially anisotropic critical behaviour (n > ric). 
Recently, H. Kleinert and V. Schulte-Frohlinde have found the RG functions for the 



(4 — e)-dimensional hypercubic model in the five-loop approximation |T^. To obtain the 



RG series of unprecedented length for the model with two quartic coupling constants these 



authors have employed the early results of five-loop RG calculations for 0(n)-symmetric y?^ 



field theory [13|. It is well known that in the replica limit n — > the scalar ip^ field theory 
with 0(r;,)-symmetric and hypercubic self-interactions describes the critical behaviour of 
the RIM provided the coupling constants have proper signs Hence, the RG expansions 



obtained in |jT2| may be used for calculation of RIM critical exponents as power series in ^/e 
or, more precisely, for extension of known three-loop expansions [^,0 up to five-loop (e^) 
terms. To find such five-loop expansions is the main goal of this paper. 

Another goal is to get numerical estimate for ric starting from the e-expansion for this 



quantity obtained in Ref. [|l^] and to compare this estimate with its analogs found earlier 
from RG expansions in three dimensions If the results given by these two approaches 

are in accord we shall be able, at last, to answer the old question: is a cubic crystal effectively 
"isotropic" at the critical point? 



II. RG FUNCTIONS AND CRITICAL EXPONENTS FOR RANDOM ISING 

MODEL 



We begin with the standard Landau Hamiltonian for a model with hypercubic anisotropy, 
describing numerous magnetic and structural phase transitions in solids. It reads: 
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where cp is an n-component order parameter, ~ r, with r = ~ " being the reduced 
deviation from the mean-field transition temperature. In the replica limit, the Hamiltonian 
Eq. is known to describe the RIM provided f o > and tto < 0. 



If we set n to zero in formulas obtained in ||T2[, we arrive, directly, at five-loop /5-functions 
and critical exponents for the model under consideration: 
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-«-^(-^-fc(3)-12C(4))-.^(-^-^C(3)-3C(4)) , (2.5) 

where numbers C(3) = 1.202056903; C(4) = 1.082323234; C(5) = 1.036927755; C(6) = 
1.017343062; ^(7) = 1.008349277 are the values of the Riemann ^-function. The coordinates 
of the random infrared-free fixed point being zeros of the above /^-functions, they are found 



in the following form |15| , p!6| : 



u 



-A^i + + Cv^' + De' + i^v^' + . . . , (2.6) 



v* = U^e+{F-p)e+iG-^C)y^' + {H-^D)e' + iL-^K)^' + ... . (2.7) 

Substituting u* and v* into the beta-functions, we get algebraic equations which decouple 
into 4 pairs of equations for A and F, B and G, C and H, and, at last, for D and L, respectively. 
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Note, that the coefficient K may be computed only in the sixth-loop approximation. After 

straightforward but cumbersome calculations one is led to the following expressions: 
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The impure fixed point has been proved to be stable in the framework of the perturbation 
theory in ^/e. Hence, we have to insert u* and v* into the critical exponents series. The 
expansions for critical exponents in powers of ^/e being the eventual goal of this study are 
as follows: 
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The A/e and e terms in Eq. (|2.10|) and first three terms in Eq. ( |2.11|) coincide with those 



calculated earlier |]15|,|16i while the rest are essentially new. Estimating coefficients in the 
above expansions numerically we obtain: 

rj = -0.0094339622e + 0.034943501v^^ - 0.044864982e2 + 0.021573216v^^ , (2.12) 

1 ■? 
= - + 0.084115823x/e - 0.016632032e + 0.047753505v^ + 0.27258431e^ . (2.13) 

Corresponding -y/e-expansion for the susceptibility exponent reads: 

7 = 1 + 0.16823164v/e - 0.028547082e + 0.078828812v^^ + 0.56450485e2 . (2.14) 

The most striking feature of the series just obtained is the quite irregular behaviour of 
their coefficients. It may be considered as one of the main reasons of failure of our attempts 
to apply various resummation techniques to these expansions. Indeed, the numerical results 
found by means of several methods based on the Borel transformation turned out both 
to contradict to exact inequalities and to differ markedly from estimates given by 3D RG 
analysis [^J^J^ and by computer simulations [17|. For example, numerical estimates for the 



exponent u thus obtained lead, via scaling relation a = 2 — Du, to positive (and big!) values 
of the exponent a, in obvious contradiction with the inequality a < proven for impure 



systems [18 



The "odd" behaviour of coefficients of the -y/e-expansions for critical exponents inherent 
in the RIM appears to be the rule rather than the exception, this apparently indicates on 
the Borel non-summability of perturbative series. The point is that the analysis of the 
perturbative expansions for the free energy of the zero-dimensional ("toy") model with 
quenched disorder has lead to the conjecture that for disordered systems perturbative series 



are Borel non-summable |19]. Such an non-summability has then been related to Griffiths 



singularities pOf . So, the irregularity of signs and values of coefficients of the -y/e-expansions 



found may be regarded as a manifestation of the Borel non-summability of these series. 



III. IS A CUBIC CRYSTAL "ISOTROPIC" AT THE CRITICAL POINT? 



Much better situation takes place in the case of a pure system with cubic anisotropy. 
Five-loop e-expansion for the marginal spin dimensionality found by Kleinert and Schulte- 



Frohlinde 12 
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= 4 - 2e + 2.58847559e2 - 5.87431189e^ + 16.8270390e^ (3.1) 

is seen to be alternating. Moreover, coefficients modulo of its Borel transform are easily 
shown to monotonically decrease what may be thought of as a manifestation of the Borel 
summability of the original series. 

Let us calculate ric for the three-dimensional system applying the Pade-Borel resumma- 
tion technique to the expansion Eq. (|3.1| ) and then putting e equal to unity. In so doing, 
however, it is very important to trace how sensitive is the numerical estimate for ric thus 
obtained to the perturbative order employed. That is why we calculate here ric not only in 
the five-loop approximation but also in three- and four-loop orders in e. Correspondingly, 
Fade approximants [1/1], [2/1] and [3/1] are used for analytical continuation of the Borel 
transforms of the original series. The results obtained in three subsequent approximations 
mentioned are as follows: 

n^^) = 3.004, n^^) = 2.918, n^^^ = 2.855. (3.2) 

These estimates are seen to behave quite regularly. They decrease with increasing the 
order of the approximation demonstrating the tendency to go deeper and deeper below 3. 
This fact is in agreement with the conclusion about the character of critical behaviour of 
cubic crystals made earlier on the base of higher-order RG calculations in three dimensions 
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||7|j9|JT0[] . Indeed, two-loop, three-loop [1,0] and four-loop |^ RG expansions for 3D cubic 
model resummed by means of the generalized Pade-Borel (Chisholm-Borel) method lead to 
the estimates: 



n 



= 3.01, = 2.95, n^^^ = 2.90. (3.3) 



These values being very close to their counterparts resulting from the e-expansion 
Eq. (|3.1|) also go down when the order of the approximation grows up, and the most accurate 
3D estimate available n^^ = 2.90 is appreciably smaller than 3. 

Pretty good agreement between the higher-order RG estimates for obtained in three 
and (4 — e) dimensions enable us to believe that both estimates are close enough to the 
exact value of n^. Hence, we may conclude that ric < 3 and cubic crystals with vector order 
parameter (3D spins) should demonstrate, in principle, anisotropic critical behaviour with 
critical exponents differing from those of the 3D Heisenberg model. On the other hand, since 
the difference between and 3 is rather small the cubic fixed point of the RG equations 
should be located in close vicinity of the Heisenberg fixed point at the flow diagram. As a 
result, critical exponents describing the anisotropic critical behaviour should be numerically 
close to the critical exponents of the Heisenberg model. 



IV. SUMMARY 

In the paper, A/e-expansions for critical exponents of the weakly-disordered Ising model 
are calculated up to the five-loop order. Coefficients of the expansions obtained are found 
to exhibit rather irregular behaviour preventing these series from to be resummed by means 
of the procedures based on the Borel transformation. This fact may be thought of as a 
reflection of Borel non-summability of such expansions conjectured earlier on the base of 
studying of much simpler (zero-dimensional) model. The marginal spin dimensionality ric 
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for a cubic model is calculated by the Pade-Borel resummation of corresponding five-loop 
e-expansion obtained by Kleinert and Schulte-Frohlinde. The estimate ric = 2.855 thus 
obtained is found to agree well with its analog extracted earlier from four-loop 3D RG 
expansions. The conclusion is made that the exact value of Uc is smaller than 3 and cubic 
crystals with vector order parameters should demonstrate, in principle, anisotropic critical 
behaviour. 
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